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Dynamical barriers of pure and random ferromagnetic Ising models 

on fractal lattices 

Cecile Monthus and Thomas Garel 
Institut de Physique Theorique, CNRS and CEA Saclay, 91191 Gif-sur-Yvette, France 

We consider the stochastic dynamics of the pure and random ferromagnetic Ising model on the 

hierarchical diamond lattice of branching ratio K with fractal dimension df = (ln(2A'))/ln2. We 

adapt the Real Space Renormalization procedure introduced in our previous work [C. Monthus 

and T. Garel, J. Stat. Mech. P02037 (2013)] to study the equilibrium time teq{L) as a function 

CO , of the system size L near zero-temperature. For the pure Ising model, we obtain the behavior 

teq{L) ^ L^e^^''^ " where dg = df — 1 is the interface dimension, and we compute the prefactor 

^^ I exponent a. For the random ferromagnetic Ising model, we derive the renormalization rules for 

^^ , dynamical barriers Beq{L) = (Inteq/P) near zero temperature. For the fractal dimension df — 2, 

5-H ' we obtain that the dynamical barrier scales as Beq{L) = cL + L^'^u where m is a Gaussian random 

*^| variable of non-zero-mean. While the non-random term scaling as L corresponds to the energy-cost 

of the creation of a system-size domain-wall, the fluctuation part scaling as L^'^ characterizes the 

barriers for the motion of the system-size domain- wall after its creation. This scaling corresponds 

to the dynamical exponent ^l> — 1/2, in agreement with the conjecture ip = da/2 proposed in [C. 

Monthus and T. Garel, J. Phys. A 41, 115002 (2008)]. In particular, it is clearly different from the 

P^ ' droplet exponent ~ 0.299 involved in the statics of the random ferromagnet on the same lattice. 

I. INTRODUCTION 

j^ I Among real-space renormalization procedures [l|, Migdal-KadanofF block renormalizations [2| play a special role 

Cj . because they can be considered in two ways, either as approximate renormalization procedures on hypercubic lattices, 
I ' or as exact renormalization procedures on certain hierarchical lattices [3|-|5|. Besides the study of pure models, these 
'^ , hierarchical lattices have been also much used to study the equilibrium of disordered classical spin models, such as 
C ■ the diluted Ising model [0] , the random bond Potts model [T-IO] , the random field Ising model [ll|, [14] and spin- 
glasses [13] • For spin-glasses, these hierarchical lattices have been also used to study dynamical properties |14| - KL7| . 
The equilibrium properties of disordered polymer models have also been considered, in particular the wetting on a 
I disordered substrate flS-Ku!] and the directed polymer model f2u-^9]. 
^~^ . In the present paper, we consider the pure and the random ferromagnetic Ising model on the hierarchical diamond 

J|T ' lattice in order to study its dynamical properties near zero temperature. We adapt the real space renormalization 
j!^ , procedure introduced in [30[ : using the standard mapping between the detailed-balance dynamics of classical Ising 
I ' models and some quantum Hamiltonian, we obtain the appropriate real-space renormalization rules for the quantum 
fvq Hamiltonians associated to single-spin-fiip dynamics. This approach has been used previo usly to characterize the 

■ ' dynamics of the random ferromagnetic chain [3Q] , of the pure Ising model on the Cayley tree [30| , of the random Ising 
r-^ ' model on the Cayley tree [S^l , and for the hierarchical Dyson Ising model [Slj . 

rf-\ I The paper is organized as follows. In section |IT1 we introduce the notations for the stochastic dynamics of Ising 

T-H ' models defined on the hierarchical diamond model. In section IIIIl we analyze via real-space renormalization the 
IL^ . dynamics of the pure Ising model. In section IIVI we study via real-space renormalization the dynamics of the 
random ferromagnetic Ising model. Our conclusions are summarized in section |Vl In Appendix \^ we derive some 
renormalization formula that are used in the text. 



o 



X 



II. MODEL AND NOTATIONS 

A. Hierarchical diamond lattice of branching ratio K 

As shown on Fig. [U the hierarchical diamond lattice of branching ratio K is constructed recursively from a single 
link called generation n = : generation n = 1 consists of K branches, each branch containing 2 bonds in series ; 
generation n = 2 is obtained by applying the same transformation to each bond of the generation n = 1. At generation 
n, the length L„ between the two extreme sites is 

Ln = 2" (1) 





n = ^ n=\ *- n = 2 

FIG. 1: Hierarchical construction of the diamond lattice of branching ratio K 



whereas the total number _B„ of bonds is 



so that 



S„ = (2X)" = Lt' 
ln{2K) 



V 



In 2 



(2) 



(3) 



represents the fractal dimension. 

As recalled in the introduction, the equilibrium of many statistical physics models have been studied on this lattice. 
Here we consider the pure and the random ferromagnetic Ising model with the classical energy 



U{C) = -J2J^JS^SJ 



i<j 



to study the properties of stochastic dynamics satisfying detailed balance. 



(4) 



B. Dynamics satisfying detailed balance 



The stochastic dynamics is defined by the master equation 

dPt (C) ^ 

c 



dt 



- ^ Pt {C) W {C -^C)-Pt {C) Wont {C) 



(5) 



that describes the time evolution of the probability Pt{C) to be in configuration C at time t. The notation W (C -^ C) 
represents the transition rate per unit time from configuration C" to C, and 



Woutic) = J2wic^c') 



(6) 



c 



represents the total exit rate out of configuration C. 

The convergence towards Boltzmann equilibrium at temperature T = ^ in any finite system 



where Z is the partition function 



Pe,iC) 



^ = E' 



-PU{C) 



.-PU(C) 



(7) 



(8) 



can be ensured by imposing the detailed balance property 

e-P^^c)W {C ^ C) - e-^^(^')W^ {C -^ C) 
It is thus convenient to parametrize the transition rates as 

W{C^C')^G (C, C) e-4[^(^')-t/(c)] 
where 



G [C, C) = G (C, G) ^ ^W (C -> C") W (C ^ G) 



(9) 
(10) 
(11) 



is a symmetric positive function of the two configurations. Near zero temperature, it is convenient to introduce the 
notion of dynamical barrier B defined by the asymptotic behavior 



G{G,G') oc e-^«('^-'^') 

/3->-+oo 



(12) 



C. Associated quantum Hamiltonian 



The standard similarity transformation (see for instance the textbooks |33l - l35l | or the works concerning spin models 



Pt{C) 



\U{C) 



^t(C) = e-^^(^)<C|V't> 



(13) 



transforms the master equation of Eq. [5] into the imaginary-time Schrodinger equation for the ket |^t > 



-\iPt>=-H\i;t> 



with the quantum Hamiltonian 



H=Y, G{G,G') [e-^N^')-^"^]|C >< G\ - \G' >< G\ 

C,C' 



(14) 



(15) 



The groundstate energy is Eq = 0, and the corresponding eigenvector corresponding to the Boltzmann equilibrium 
reads 



e-^uiC) 



l*o>^E^7^|c> 



(16) 



The other energies En > determine the relaxation towards equilibrium. In particular, the lowest non-vanishing 
energy Ei determines the largest relaxation time (l/Ei) of the system, i.e. the 'equilibrium time' needed to converge 
towards equilibrium, 



_ 1 



(17) 



D. Single-spin flip dynamics of Ising models 



We will focus here on single spin-flip dynamics satisfying detailed balance of Eq. [9l where the transition rate 
corresponding to the flip of a single spin Sk reads 



W (Sk ^ ~Sk) ^ G"" 



hk — 2_^ JikS'i 



,-ps4j:,^,j,,s,] 



(18) 



G™'[/ife] is an arbitrary positive even function of the local field hk = J2id^k Jik^i- For instance, the Glauber dynamics 
corresponds to the choice 



^Glauber W 



2cosh(/3/i) 



(19) 



The quantum Hamiltonian of Eq. [15] reads in terms of Pauli matrices (cr^, cr^) [30|, I36l44: 



n =2^G' 






-/3^^(E.^fcJ.fc'^0 



<) (20) 



For ferromagnetic models near zero-temperature, more precisely when the temperature is much smaller than any 
ferromagnetic coupling Jy 

< T < Jy (21) 

the thermal equilibrium is dominated by the two ferromagnetic groundstates where all spins take the same value, and 
the largest relaxation time t^q — l/Ei corresponds to the time needed to go from one groundstate (where all spins 
take the value +1) to the opposite groundstate (where all spins take the value —1). The aim of the renormalization 
procedure introduced in [30| is to preserve the lowest non- vanishing energy Ei of the quantum Hamiltonian. We 
have already explained the application of this renormalization to the random ferromagnetic chain (30l] , to the random 
ferromagnetic Cayley tree [32], and to the Dyson hierarchical Ising model jSliJ. In the following, we derive the 
appropriate renormalization rules for the diamond hierarchical lattice described above. 

III. DYNAMICS OF THE PURE ISING MODEL 
A. Principle of the Renormalization for the dynamics 

To analyze the statics of the ferromagnetic Ising model on the diamond lattice, one has to renormalize bonds to 
obtain the RG rules for renormalized couplings. However here to analyze the single-spin-flip dynamics, we wish to 
define renormalized spins that represent ferromagnetic clusters of spins flipping together. 

We start from the diamond lattice with n > 2 generations, of length L„ = 2", containing i?„ = {2K)"' bonds. All 
bonds have the same initial ferromagnetic coupling J. All elementary spins have a dynamics characterized by the 
function G'^^'lh] of the local field h {Eq^. 

After one RG step, we wish to have a diamond graph of generation (ti — 1), where there is an alternation of initial 
spins, with a dynamics still characterized by the function G™*[/i], and of renormalized spins Sr-^ (representing clusters 
of {2K + 1) initial spins), whose dynamics is represented by some renormalized amplitude Gfn. The renormalized 
bonds have for renormalized couplings Jr-^ — KJ . 

More generally, after p RG steps with 1 < p < n— 1, we wish to have a diamond graph of generation (n — p), where 
there is an alternation of initial spins, with a dynamics still characterized by the function G""*[/i], and of renormalized 
spins Sr (representing clusters of {2K + 1)^ initial spins), whose dynamics is represented by some renormalized 
amplitude Or . The renormalized bonds have for renormalized couplings 

Jr, = KPJ (22) 

In particular, after p ~ n — 1 RG steps, we wish to have a diamond graph of generation 1, where the two bound- 
ary spins are initial spins, with a dynamics still characterized by the function G*"*[/i], and where the K internal 
renormalized spins Sr^_-^ (representing clusters of {2K + 1)"~^ initial spins), whose dynamics is represented by some 
renormalized amplitude Gr^_j^. The renormalized bonds have for renormalized couplings Jr^_i = K"~^J. 

The last RG p — n RG step is thus special : the {K + 2) remaining spins are grouped together into a single renor- 
malized spin Sr^ (representing the whole sample of {K{2K + 1)"~^ -I- 2) initial spins) whose dynamics is represented 
by some renormalized amplitude G^"*. 

Let us now explain how the RG rule for the renormalized amplitude Gr describing the dynamics can be obtained 
for the special last RG step p = n and for the bulk RG steps 1 < p < n — 1 respectively. 

B. RG rule for the last RG step p = n 

The last RG step p = n shown on Fig. [2] involves two boundary spins Sei and S'e^ and K internal spins 
(•Sai , •S'aa, -SaK) with the foUowing classical energy (Eq. U]) 

K 
U{C) = -Jr,^_, Y, SaASe, + Se,) (23) 




Si 



ast 



FIG. 2: Last RG step for the dynamics of the pure Ising model : notations for the quantum Hamiltonian of Eq. [24] 
The quantum Hamiltonian of Eq. [501 associated to the single-spin- flip dynamics reads 

K 



H 



K+2 



G' 



+G'' 



K 






K 



-P-jRr,-i<S''l^+<2'> 



<) 



We are interested in the lowest non-vanishing eigenvalue Ei > 0. The eigenvalue equation 

^ iHK+2 - Ei)Ui > 
for the corresponding eigenvector in the basis 

\lpl>= ^ ^ X! ■• X! '^l{Sei,Se2,Sai:SaK)\Sei,Se2,Sai:SaK > 



reads 



K 






=[G"" Jfl„_,5^5a, 
i=l 
K 

+G™* Jr„_,J2^-^ 

i=l 
K 

^l( i>ei 1 '-'62 1 ^ai ■■^aK/ 
rlv'-'ei 1 ^e2^ ^ai ■■^aK) 



-G* 



-C 



K 



' 


K 


Jr^ 


-.>:s". 


_ 


i=l 


- 


if 


Jr^ 


-.>>. 


. 


i=l 



(24) 



(25) 



(26) 



(27) 



i=l 



Let us now use the symmetry between the K spins Sai to note ^(-S'ei , k, Se^ ) the components i/)i {Se^ , 5*62 , Sa^.-Sa^ ) 
where k € 0, 1, .., ii' spins among (S'^j , .., Sa^) take the value (— ) 



K 



'^ ^ 1 — S* 

i=l 

Eg [27l becomes 

= [G™ [Jfl„_i(i^ - 2A:)] e-/3JH„-iS.,(K-2fc) _^ ^^^^ [jfl„_,(X - 2fc)] g-'^'^^^-i'^^sC^-sfc) 

[K - k)GR^_, [jR.^_ASe, + ^ej] e-^^«"-i(^=i+^=^) - £;i]0i(5,,, fc, 5eJ 

-kGR^^_^ [j_R,„_i(5'ei +'S'e2)] 4'liSei,k- IjS'e^) 

-(i^-fc)Gfl„_, [Jfl,„_,(5ei+5e2)]</'l(5ei,fc + 1,5^2) (29) 

For _Eo = 0, the (non-normahzed) groundstatc is known to be given by Eq. [12] with the classical energy of Eq. [531 

MSe„k,Se,) = e#-^«„-i(^-2fe)(s=i+Se2) (30) 

so it is convenient to look for the solution of Eq. [27] via the amplitude 

that satisfies 

= [G™ [Jr^_AK - 2fc)] e-''^«^-i(^-2fc)s., ^ (^«m [Jr^^AK - 2fc)] e"'''^"--! (-^-2^)^=2 

+ (i^- fc)Gfl„_, [Jfl,„_i(5ei +5^2)] e-'3''«"-i(^=i+^=^) -^l]A(5ei,fc,5e2) 
-G™ [ Jfl,„_, (i^ - 2fc)] e-'P-Jr^r.-i(K^2k)S._, A{-S^, , k, Se, ) 

-G™ [Ji^„_,(i^-2fc)] e-^^«"-i(-^-^'=)^'=2A(S'ei,fc,-5e2) 

-(i^-fc)Gfl„_, [Ji?„_i(5ei +5,2)] e^''^«"-i(^=i+^^^U(5e,,fc + l,5e2) (32) 

Let us now consider the dynamical path associated to the swap of a domain-wall between the two ferromagnetic 
groundstates with the following notations for configurations 

Co ={Se,=+l,0,Se,=+l} 

Cl ^{Se,^-l,0,Se,=+l} 

C2 ={5ei=-l,l,5e2=+l} 

Cp = {S'ei = -l,p- 1,5*62 = +1} 

Ck ={Se,=-l,K-l,Se,=+l} 

Ck+1 = {"Sei = -l,iir, 5e2 = +1} 

Ck+2 ={5e, =-l,i^,5e2=-l} (33) 

The physical meaning is that the transition from Go to Gi corresponds to the entrance of a domain- wall at the boundary 
ei, the transitions between the configurations (Gi, G2, .., G^+i) correspond to the displacement of this domain-wall, 
and finally the transition from Ck+i to Ck+2 corresponds to the exit of the domain- wall at the boundary 62- 
The corresponding amplitudes A{Cq) for this dynamical path satisfy (Eq[32|) for 2 < q < K 

= {{q- l)Gfl„_, [0] + {K + l- q)GR„_, [0] - ^1) A{Cg) 

-{q - 1)Gr^_, [0] A(G,_i) -{K + l- q)GR^_, [0] A(G,+i) (34) 



and ioTq~Q,q = l,q = K + l and q = K + 2 

= [G^™ [KJr,.-^] e-'3^'^«"-i - Ei]A{Co) - G™ [ii:Jfl„_J e-''^^''--iyl(Gi) 
= [G'"' [XJfl„_ J e^^^«"-i + KGr^_, [0] - i?i]A(Gi) 

-G"^ [if Jfl„_J e^^-^^-MlGo) - ifGfi„_, [0] A{C2) 
= [G"'' [XJfl„„ J e^^^«"-i + A'Gfl„_, [0] - E,]A{Ck+i) 

-G™ [if Jfl„_J e^^^--.-iA(GK+2) - KGr^_, [0]A{Ck) 



= [G™ [if Jfl„_ J e- 



-/3KJfl 



Si]A(Ga-+2) - G"" [KJr„_,] 



-KPJb 



.-M(Gk+i) 



(35) 



For this effective one-dimensional problem with the effective transition rates 

Wf^iCo -^ Gi) = G"" [if Jfl„_J e-^-^-'^-i 



W^^^^(G, ^G,+i) =(if + f-g)Gfl„_, [0] for 1 < g < if 



PkJr^ 



and 



W^ff{CK+i ^ Ck+2) = G™ [if Ji^„_J e 



W'=^^(Ci ^ Go) = G'"' [if Jfl„_J e'^-^'^"-! 
W^"^^(G, ^G,_i) = (g - f )Gfl„_, [0] for 2<g<if + l 
W^ff{CK+2 ^ Gk+i) = G'"' [if Jfl„_J e-^'^-^«-i 



(36) 



(37) 



we may use Eq. IA21I of the Appendix to obtain the renormalized amplitude G^^f^^^ = G]i{Cq,Ck+2) when the 
domain-wall enters by the boundary ei and exits by the boundary 62 



62 



|[(7(Co)-C/(Cic+2)] 



G 



(ei,e2) 
Rjilast 



WffiCo -> Gi) 



G'- [KJn,,_,] ^ Gr„_, [0] 



rn— 1 g- 
^2/3KJh„ 



.liM^e//(C^^C^^^. 



if 

E 



(to- l)!(if -to)! 



(38) 



Taking into account the other case where the domain- wall enters by the boundary 62 and exits by the boundary Ci, 
which actually gives the same contribution 



-<(ei,e2) 



~((e2,ei) 






1(61,62) 



i(e2,ei) 



-.(61,62) 



we obtain that the final total amplitude G^f = G'^l\2'^ + G'£\l\'^ = 2G'^l\l^J^ reads 



.P^jRr. 



r^last 

Rn 



,^pKjR„-.^^k\{K-i-ky. 



G™ [if Jii„_J 2Gr^_, [0] ^ if! 



(39) 



(40) 



C. RG rule for the bulk RG steps 1 < p < n - 1 

The bulk RG step p shown on Fig. [3] involves two external spins S*^^* and S'g^*, and (2if + 1) internal spins 
(Sai , 5*02 , ■, Saii , Sc, Sbi ,Sb^,., Sbii ) with the following classical energy (Eq. g]) 



K 



K 



u{c) = ~Jr^_, Y. Sa^ (s, + ST) - Jfl,-i E ^b. [Sc + ST) 



(41) 




/fi 



»Sr 



FIG. 3: Bulk RG step for the dynamics of the pure Ising model : notations for the quantum Hamiltonian of Eq. [42] 



The quantum Hamiltonian of Eq. [20] associated to the single-spin- flip dynamics reads 

K 



H. 



2K+1 



G* 



K 



Jfl,_i^«+0 



(' 



-/3JHp_i<T-Ef=lK.+'^6 






(42) 



Let us now focus on the external Domain- Wall conditions 



oext 



cext 



(43) 



and take into account the symmetry between the K spins Sa^ , and the symmetry between the K spins Sb^ to 
note 4>{ka,ki,) the components of ?Ai where ka G 0,1,.., K spins among {Sa^, .., Sa^) take the value (— ), and where 
ka € 0, 1, .., K spins among {Sai , .., Sa^) take the value (— ) 

^ 1 - S* ^ 1 - S* 

MST' = -1, Sa„.., Sa^,Sc, Sb„.., Sb^,sT = +1) = Mka = Y. -^^ s„ kb^Y. -^) (44) 



Then the eigenvalue equation = {H2K+1 — £'i)|V-'i > reads 

= (G™ [Jb.,_,{2K - 2ka - 2kb)] e-'3''Hp-i(2-ft^-2fco-2fc,)S. 
+kaGR^_, [Jr^_,{S, - 1)] e^'^«.-(^-i) 
+ iK - K)Gr^_, [Jr^_ASc - 1)] e^^-'H.-i^^-i) 
+fc6G^,_.[JH,_.(5. + l)]e^-^-.-(^^+i) 

+ {K - h)GR^_, [Jr^_, {Sc + 1)] e-''''-p-i(^=+'^ - Ei)Mka,Sc, h) 
-G™^ [Jr^_,{2K - 2ka - 2h)] Mka,-Sc,kb) 
~kaGR^_^ [Jr^_ASc - 1)] (t>iika - l,Sc,kb) 
^{K - ka)GR^_, [Jr^_, {Sc - 1)] Mka + 1, 5c, U) 

~kbGR^_^ [Jrp_JSc + 1)] <i)i(ka,Sc,kb - 1) 

-{K-h)GR^_, [jR^_ASc + l)\4>i{ka.ScM + l) (45) 

For _Eo = 0, the (non-normahzed) groundstate is known to given by Eq. [16] with the classical energy of Eq. [41] 

Mka,Sc,kb) = el^«P-l[Sc(2K-2fc„-2fc,) + (2fc„-2fc.)] ^^g^ 

so it is convenient to look for the solution of Eq. [27] via the amplitude 

A(fca,5c,fc.) = ^4H^ (47) 

that satisfies 

= (G™' [Jr^_,{2K - 2A:a - 2kb)\ e-0-'i^,-i(2K-2k,-2k,}S^ 

+kaGR^_, [Jr^_ASc - 1)] e^^-p-(^-i) + (i^ - fc.)G;,,^_, [J«,_,(5c - 1)] e-^^«.-(^-i) 

+fc;,Gfl^„, [Jfl,,_i(5c + 1)] e'3'^«P-i(^=+i) + (i^- fcb)Gfl^_, [Jr,-ASc + 1)] g-'^^^p-i^^^+i) - Ei)A{ka,Sc,kb) 

-G™* [Jfl^_,(2i^- 2fc, - 2fc,)] e-^^'^p-i(2^-2'=»-2'=')^=A(A;„-5e,fcb) 

-fcaGi,^_, [Jfl,_i(5c - 1)] e^-'---i(^^-i)A(fc, - l,5c,fc,) 

-(/f - ka)GR^_, [Jr^_ASc - 1)] e~'''''^.-i(^^~^)A(fc, + l,Sc,kb) 

-kbGR^_, [Jr^_ASc + 1)] e''-''''"^'-^^+'^A{ka,Sc,kb - 1) 

-(i^-fcb)Gfl^_, [jR^_,{Sc + l)]e~^'''''-^'^^^+''^A{ka,Sc,kb + l) (48) 

Let us consider the dynamical path along the following (2i4r + 2) configurations 



Go 


- {0, 5c = 


= +1,0} 


Gi 


= {l,5c = 


= +1,0} 


G2 


= {2, 5c = 


= +1,0} 


Gk 


= {i^, 5c = 


-+1,0} 


G/<-+i 


= {K, Sc 


= -1,0} 


G/i:+2 


= {K, Sc 


= -1,1} 


Ga:+3 


= {K, Sc = 


= -1,2} 


C2K 


= {K, Sc 


= -l,i^-l} 


C2K+1 


= {K, Sc = 


= -l,K} 



(49) 



in order to describe the motion of a domain-wall corresponding to the boundary conditions of Eq. 
Then Eq. [45] becomes for < q < i^T - 1 



= {qGR^_, [0] + {K- q)GR^_, [0] - i?i)A(G,) - qGR^_, [0] A(G,_i) - {K - q)GR^_, [0] A(G,+i) (50) 
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foiq — K 

= (G™ [0] + KGb,_, [0] - Ei)A{Ck) - G™ [0] A(Gk+i) - KGr^_, [0] A(Gk-i) (51) 

ioi q = K + 1 

- (G™ [0] + KGr^_, [0] - £;i)A(Gk+i) - G™ [0] A(Gk) - KGr^_, [0] A(G;f+2) (52) 

and ior q = K +l + k ioi I < k < K 

= fcGfl,^_, [0] + {K- k)GR^_, [0] - ^i)A(GK+i+fc) - kGR^_, [0] A(GK+fe) - (i^ - k)GR^_, [0] A(GK+2+fc)(53) 
For this effective one-dimensional problem with the effective transition rates 

W^^-^-^(G, -> G,+i) ={K-q)GR^_,[Q] for < g < i^ - 1 



W^^iCx -^ Gk+i) = G'"' [0] 
W^^-^^(G, ^ G,+i) = (2/^ + 1 - q)Gi^^_, [0] for K + l<q<2K 



and 



W^ff{Gq^C^-i) =qGR^_,[Q] for l<q<K 

W^^iCK+i -^Ck) = G™* [0] 
W^ff{Gq^Cq-i) ={q-K-l)GR^_,[(S\ for K + 2<q<2K + l 

we may use Eq. IA21l of the Appendix to obtain the renormalized amplitude Gr [0] = Gr{Gq, G2K+1) 



Gr„ 



g|[C/(Co)-t/(C2K + l)] 

VF«//(Go ^ Gi) 
1 2 



^"^ (m)!(i^-l-m)! 



G""[0] ^Gfl^_,[0]^^^ (i^-1)! 



E 



(54) 



(55) 



(56) 



D. Conclusion 

Using the combinatorial formula [43] concerning the sum of the inverse of binomial coefficients 

^"^ 'm)\{K -l-m)\ ^^ 1 _ _^ ^ 2'= 



E- 

m=0 



{K-l)\ 



"-1 /c=l 



and the form of renormalized couplings of Eq. 



./i?„_, = i^"" V 



we may rewrite the last RG step of Eq. |40] as 



1 



aPK^J 



„2/3_fS'"J K ^f, 



V- 

Qlast Q^n^ [^nj] 2^+lGfl„_, [0] f^ fc 



+ 



fc=l 



and the bulk RG steps 1 < p < ?i — 1 of Eq. [55] as 

1 1 



2fe 



Gfl, [0] G™^ [0] ' 2^'-iGfl^_, [0] ^^ fc 



with the initial condition 



(57) 



(58) 



(59) 



(60) 



Gr, [0] = G"" [0] 



(61) 
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In terms of the numerical constant 



^^-:^Et (62) 



2K-1 A^ k 

the solution of the recurrence of Eq. [^01 reads for < p < n — 1 

1 1 ^ 

G^^=G^^E(^-)' (63) 

so that the final renormalized amplitude of Eq. [59] for the whole sample reads 

1 pl3K"J 2I3K'\J « 

+ I?5i^E(--)' (64) 



^tn L J L J (7—1 



The final conclusion is thus that the equilibrium time tj^ needed to go from one groundstate (where all spins 
take the value +1) to the opposite groundstate (where all spins take the value —1) reads for the Ising model on the 
diamond lattice of branching ratio K with n generations reads 

1 pl3K"J „2l3K^J " 

'" ~ G'^^* G™' [K^J] 4G™i [0] ^^''^> ^"^^ 

In particular, the dynamical barrier B^"' defined by the low-temperature exponential behavior reads in terms of the 
length L^ and of the fractal dimension dj (Eqs [U and [3] ) 

S(«)= hm ^ = 2JX" = 2JLf/"' (66) 

in agreement with the expected scaling U^' of the energy cost of an interface of dimension ds = d — 1 \n a, space of 
dimension d. 

Here, we have in addition computed explicitly the prefactors in Eq. [55] : the leading behavior is a power law of the 
length Ln 

df-l df-1 

"eg 

where the exponent reads in terms of the constant ck of Eq. 



4;) ^g/32JL/- (^^)«^g/J2JL/- ^OK (g7^ 



OiK 



\ncK ^"^ 



1 ^K 2i 

2^-1 Z^fc=l fc 



(68) 



In 2 In 2 

In particular, for the first values of K, we obtain 

aK=i = 1 
aK=2 = 1 

In I 
"^=^ ^h^ 

In I 

For K = 1 corresponding to the one-dimensional chain, the result aK=i = 1 is in agreement with previous studies 
(see [201 and references therein). 

IV. DYNAMICS OF THE RANDOM FERROMAGNETIC MODEL 

A. RG rule for the last RG step p — n 

The last RG step p — n shown on Fig. [3| involves two boundary spins Se^ and Se2 and K internal spins 
{Sai, Sa2, ••S'aic) with the following classical energy (Eq. [4]) 

K 

U{C)^~J2 ^-^ (-^-^ ^-1 + J'a, ^-- ) C^O) 

i=l 
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Si 



ast 



FIG. 4: Last RG step for the dynamics of the disordered Ising model : notations for the quantum Hamiltonian of Eq. [71] 



The quantum Hamiltonian of Eq. [20] associated to the single-spin- flip dynamics reads 
Hk+2 = G'"' 



■ K 



+G'' 



.4=1 

K 



K 

E 

i=l 



i=l J 



(71) 



Let us first consider the dynamical path where the spins are flipped in the order Se^, Sai, Sa2, ■■, Sa^j Se2- It is 
convenient to introduce the following notations for the corresponding {K + 3) configurations 

Go = i'-'ei = +li "Jai = +1; >-'a2 = +1; ••j ^ax ^ +1? '^£2 = +1/ 

Ci = jOei = ^Ij "Jai = +1; '-'02 = +1) ••! "-"qk' ^ +I7 '-'62 ^ +1/ 

C2 = {Sei = —1, 5*01 = —1, Sa2 = +1, •■, Saii = +1, >S'e2 = +1} 

^P ^ i'-'ei ^ ^-^1 "^ai ^ ^1; ■■■! '-'ap-i ^^ l; '-'op ^ +1.., "Jafc ^ T^-I^j '-'e2 ^ ' ^i 

L-Zf = i'-'ei = ~li "JQi — ^1; ■•■: "-'a/c-i = ~1: j '-'ttK ~ ^^^l '-'62 ~ ~^"1| 

Ck+1 ~ {Sei = ^1, Sai ~ —1, ..., "S'a^ = —1, Se2 = +1} 

Gk+2 = {Sei = ^1, "Sqi = ^1, ■•■, '5'a/c — ~lj 'S'e2 = ^1} (72) 

The physical meaning is that the transition from Go to Ci corresponds to the entrance of a domain- wall at the boundary 
ei, the transitions between the configurations (Ci, C2, .., C/f+i) correspond to the displacement of this domain- wall, 
and finally the transition from Ca'+i to Gk+2 corresponds to the exit of the domain-wall at the boundary 62. 
The two boundaries configurations are the two ferromagnetic groundstates of classical energy (Eq |70p 



K 



[/(Co) = UiGK+2) = -Y.(Ja. + J'a.) 
i=l 

The intermediate configurations Gp for 1 < p < K + I have a higher classical energy (Eq [70l) 



p-i 



K 



U{Gp) = J2^-Ja. + J'aJ + Y^iJa. " J'aJ 



(73) 



(74) 



So we may used the formula of Eq. I A23I derived in the Appendix for the renormalized amplitude along this dynamical 
path 



^^(ei,ai,a2,..,aif,e2) z,^ /-< \ 

'-'R \'^0,'~^K+2) 



g-|[C/(Co)+C/(CK+2)] y^ ^' 

m=0 G {G,n,G„i+l) 



(75) 
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Taking into account the local field on the spin that is flipped between two consecutive configurations, one obtains for 
the two boundary flips 





r K 1 


ini 


>>.. 




Li=l J 




[ K 1 


ini 


>y; 




Li=l J 



G(Co,Ci) =G 

G {Ck+i,Ck+2) =G 
and for the intermediate flips 1 < m < K 

G (Cm, Cm+i) = C^a„, [Ja,„ ~ Ja„A 

The renormalization formula of Eq. [7S] for the last RG step finally reads 



{ei,ai,a2,..,aK,e2) 



G 



Rlast 



JT.^=iJ' 



K 



G™ Y.t=lJ'a. 



(76) 



(77) 



(78) 



B. RG rule for the bulk RG steps 1 < p < n 1 




'R 



o5, 



J\ 



R 



5exi Q, ext 

B "^S 

FIG. 5: Bulk RG step for the dynamics of the disordered Ising model : notations for the quantum Hamiltonian of Eq. |80] 
The bulk RG step p shown on Fig. [5] involves two external spins S*^^* and 5*^^*, and {2K + 1) internal spins 



(S'ai , 5*02 , ., SaK , Sc, Sbi , 5*62 , ., Sb^ ) with the following classical energy (Eq. g]) 



K 



K 



uic) = -J2^-^iJ-^^- + JL^a') -Y.^bAJb.^c + Jfc.^lf') 



(79) 
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The quantum Hamiltonian of Eq. [20] associated to the single-spin- flip dynamics reads 

" K 



H- 



2K+1 



G 






K 



'^b. 



(80) 



Let us first consider the dynamical path where the spins are flipped in the order Sai , Sa2 , ••, Sa^^ , Sc, Sbi , 5*62 ? ■■! "5*6;, 
More precisely, for the external Domain- Wall conditions 



Qext 
next 



= -1 

= +1 



(81) 

it is convenient to introduce the following notations for the corresponding (2K + 2) configurations that describe the 
motion of the domain-wall 

Co ~ {Sa^ = +1, Sa2 = +1, •■, SaK = +1, Sc = +1, Sb^ = +1, Sfc^ = +1, .., 5*6^^. = +1} 

Ci — {Sa-^ = —1, Sa2 = +1, •■, SaK = +li Sc = +1, Sb^ = +1, Sfc^ = +1, .., 5*;,^^. = +1} 

C2 — {Sai = —1, Sa2 = ^1; Sas = +1.., Sa^ = +^, Sc — +1, Sbi = +1, Sb2 ~ +1, ••, Sbji = +1} 

\Sai = ^Ij '-'02 = ^1; ■•! Saii = ^1, Oc = +l5 'Jfci = +lj '-'62 ^ ^"l' •■> '-'bfc ^ ^1| 



C/f+1 

Ck+2 
Ck+3 

C2K 
C2K+1 



{Sai 
{Sai 
{Sai 
{Sai 



I7 '-'02 
I7 '-'02 
I7 '-'02 
I7 '-'02 



— {Sai — 17*502 — J^7 



7 SaK — —1, Sc — — 17 S'bi — +I7 Sb2 — +I7 ••! "SfcA' — +1} 

7 Saj^ = ^1, Oc ^ ^17 Sbi ^ ^^7 '-'62 ^ T"^7 •■! '-'bA' ^ ''-I 

7 '-'ajc ^^ ^li "Jc ^= I7 »-'6i ^^ ^17 '-'62 ^ ^^7 Sb3 = +1. •■7 t^bji ^ 

7 Saj^ = ^1, Oc = ^17 »J6i = ^17 '-'62 ^ ^7 •■»-'6a'_i ^ ^17 »56a- 



-1} 



-1} 



7 SaK — ^17 Sc — — 17 'S'bi — — 17 Sb2 — — 17 --SbK-i — ^^7 '5&A- — ^1} (82) 



The classical energy (Eq [7^ of these configurations read for < to < ii' 



K 



K 



U{C^ 



Y.^Ja.-J'a.)~ J2 (Ja.~J'aJ'Y.(Jb^+JL) 



i=l 
K 



K 



UiCK+l+,n) =Y.^-Ja^-J'a^)+Y.^~Jb^+Ji)- E (--^"^ +'^U (83) 

i— 1 z— 1 i— m+1 

So we may used the formula of Eq. IA23I derived in the Appendix for the renormalized amplitude along this 



dynamical path 



G(ai,a2,..,aA,c,6i,b2--,fcA)//-< /-< \ 



- „-f [C/(Co)+C/(C2A + l)] 



2^ g|[;7(C,„)-(-t/(C„+i)] 



E 

m=0 



G [Gm, Gm+l) 



(84) 



Taking into account the local field on the spin that is flipped between two consecutive configurations, one obtains for 
the two boundary flips for 1 < ^?^ < -^^^ ^ 1 

G [Cm, Gm+l) — Gam ["^a™ ~ ''a^J 



G{Gk,Gk+i) = G* 



K 



J2^Ja. - Jb, 



1=1 



G {GK+m,GK+m+l) — Gb,„ [Jfe^ — Jfe^^] 

For the renormalized spin Sr, the two renormalized ferromagnetic couplings read 

K 

Jr =2Z '^'°'i 



(85) 



J'« 



i=X 
K 

E4 



(86) 
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and the absolute value of the renormalized local field hfi in the domain- wall configurations takes the single value 

K 

\hR\ = \JR-fn\ = \Y.'^J'a^~JU\ (87) 

So one does not have to renormalize a function G[h] of an arbitrary local field h, but only a set of three correlated 
variables {Jr, J'^, Qy^^^°-'^^-^°-'<,c, i, 2--, k) j^^ _ j^ _ j^j^ representing the two ferromagnetic couplings {Jr, J^) of Eq. 
[551 and the corresponding numerical amplitude (^^I'^^v.aic.c, i, 2--, '^'•'[/j^ — j^ _ j^] given by the final formula ( Eq. 

M) 



G^ [hR ^ Jr- Jj^l „j^i ^a™ [•Jam "'a^J 



m=l '^bm [Jbm - JbJ 

C. Analysis of the random ferromagnetic chain (K = 1) 

For the random ferromagnetic chain corresponding to the branching ratio K = 1, the renormalization formula for 
the last step (Eq. [75]) becomes 

(89) 



For the bulk, we may rewrite Eq. [551 as 

eP[j'a,+Jl,] e^'I^-i+^-il e^(^»i+-'''i) e'^[^''i+^^i] 

''[■-'ai+'-'t ] 

to make clearer that the combination J r ,, _t/ i has a simple renormalization rule. By iteration, we finally obtain 

that the final amplitude Gi^lf^ {Ln = 2") for a system of size Ln ~ 2" of n generations reads in terms of the initial 
function G™'[/i] (with the notations J-i.o = = Jl„,l„+i) 



2^ (^mj [ 7 , . _ 7- - , , 1 *- -* 



G,tf'(L„ = 2") tt G-M^.-i. - Jm+i] 

when the domain-wall enters by the boundary ei and exits by the boundary 62 

We should now take into account the other case where the domain-wall enters by the boundary 62 and exits by the 
boundary ei, which actually gives the same contribution 

GtT\L„ = 2") = Gtr\Lr.-2n (92) 

The total amplitude Giast{Ln = 2") is the sum of these two contributions 

G,MLn = 2") = Gtr\L^ = 2") + Gtr\L,, = 2") = 2G|:--)(i„ ^ 2") (93) 

so that the final result reads 



9 Z^ H'^m \T_. , , _ 7, ., 1 (^ ) 



GlastiL„=2^) 2 ^ G"" [J,_i., - J,,,+ i] 

in agreement with the results of Ec^. (104) obtained in our previous work [30] via the Boundary Renormalization 
procedure. This agreement shows the validity of the Bulk Renormalization procedure within the domain- wall approx- 
imation. 
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D. Dynamical barriers for branching ratio K > 1 



For K > 1, we have to compare the various dynamical paths that display different dynamical barriers as a conse- 
quence of the disorder. 

As explained before Eq. I88[ the important quantity is the numerical amplitude Gf{[hfi — Jr — J'ji] which is correlated 
with the two renormalized ferromagnetic couplings {Jr,J'j^). It is thus convenient to introduce the corresponding 

dynamical barrier B^^ ^' " defined by the low temperature behavior 



/3— >CXD 



,9S</^"'«' 



(95) 



where the notation B^ "' "- has been chosen to remind that this barrier is correlated with the two couplings {Jr, J'r)- 
Let us now focus on the Glauber dynamics of Eq. [19] with the following low-temperature behavior 



^GlauberW 



1 



-I3\h\ 



2cOSh(/?ft,) ;3->oo 



(96) 



Optimization of the dynamical path for the last RG step p = n 



In terms of dynamical barriers, Eq. [78] with Eq. [96] yields that the final dynamical barrier Bi^^^ 
to the given dynamical path (ai, 02, .., a^) reads 



{ai,a2,...aK) 



associated 



B 



{ai,a2,..,aK} 
last 



K 



,('^a™,-'a,J 



K 



2Y,Ja-.2Y,J'ai, ^ max^ B^T" "'" V 2 ^ J^ + J,„ + J^^ + 2 ^ J,, 



2^771+1 



(97) 



We now have to consider the K\ possible dynamical paths : for a given permutation vr of the K renormalized spins, 
the dynamical barrier ;B;^J/^'' ^'■^^'"" "'^^' associated to the path (aTr(i), 07^(2)7 --i '27r(x)) reads by adapting Egj^ 



fi, 



(a,(i),..,a^(jc)) 



last 



K 



K 



2E^-;25:j;;^max^ <:t 



(■^".(™,^^",(™,) 



K 



2E'^^.(, + ^'^.<™) + '^^.(™,+2 E •^«.<.. 



i=l 



i=l 



i=l 



2—771-1-1 



(98) 



We now have to choose the dynamical path, i.e. the permutation ir leading to the smallest barrier. So the final 
renormalized barrier Blast is given by the minimum of Eq. 1981 over the K\ possible permutations 



Blast = mm ( Bi^J 



(99) 



= mm max 



K 



K 



2yj,,;2Vj^.; max ^a"!*; 

Z^ Z^ Qj l<77l<i<- \ ^ > 



(^».(™,^^:,(™,) 



K 



2E^U,+'^'^^<™)+'^^.(^)+2 E '^'^.(.. 



2—771+1 



2. Optimization of the dynamical path for the bulk RG steps 1 < p <n — 1 



Similarly for the bulk RG step, Eq. [55] yields that the renormalized dynamical barrier gt^i'"2,-,aA',c, 1, 2-, k) 
correlated with the renormalized couplings (Eq. 



K 



Jr 


i—l 


J'r 


K 



(100) 



i=l 
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and associated to the dynamical path (ai, 02, .., Qk, c, 61, 62.., &if ) reads 



Jai,a2,..,aK,c,bi,b2..,bK) 
R 



Do ' ' = max( max 



l<m<K 
K 






K 



■i^m+1 






(101) 



i=l 



max 

\<m<K 



B. 



iJi'n^-JD 



7n — l 



K 



E (^^> - ^^.) + Jb^. - JL + E (2J&. - ^:. - 4 J 



i—m+l 



We now have to consider the (Kl)'^ possible dynamical paths : for a given permutations (tTq) of the «„ 
renormalized spins, and for a given permutations (tti,) of the 6,„ renormalized spins, the dynamical barrier 



i? — — — — 

reads by adapting Eq llOll 



associated to the path {a^^^i),a^^(^2), ■■,a^^{K),c,K,{i),Kti2), -^KtiK)) 



,(a,r„(l),ax„(2),--,ax„(K):C,6„(i),fc^ (2), ..,&„. (K)) 



Bj^ " ' "' ' " ' "• ' "• ' "■ ' = max( 



max 

l<rn<K 



K 



Ba 



(''..i^yJ'^..,^,^ 



K 



E (^-^a^^w - J'a^^^,^ - JL(^)) + J<^^ai^) ~ 4,, („) + E (^a,,„(.) "^L („) 



2^771+1 



K 



Y.^Ja. - Jb.)\ + Y.^Ja. + •^fc. - 4 - 4J; 



(102) 



i=l 



max 

l<m<K 



B 






if 



E (-^l.., - Jl..) + <(^, - JL^.. + E (2^^,.., " J^,, - ^L,<.) 



z— m+l 



We now have to choose the dynamical path, i.e. the permutations (7ra,7rt,) leading to the smallest barrier. So the 

final renormalized barrier Bj^"^' is given by the minimum of Eq. 11031 over the K\ permutations (tt^) and over the 
K\ permutations {'Kb) 

B^'^'''-^ ^ min (^(--(i)-'^-(=)-^'^-(-)'^'''".(i)'''^.(^) '■■'"'.(-))) (103) 

To see more clearly the structure, let us now focus on the case K = 2. 



E. Case K = 2 corresponding to the fractal dimension df — 2 

1. Last step RG rule for dynamical barriers when K — 2 
For K = 2, Eq. 1971 involves a maximum over four terms 



B 



(01,02) 



last 



2{Ja, + Ja^y, 2{J'a, + J'aJ; Ba^'^""'' + Ja, + J'a, + ^Ja.^Ba:''" "'' + 2j'„, + J,, + j/ 



'ai 

and Eq. IIOOI involves the minimum over two permutations 
Blast = min( max 

2{Ja, + Ja.y 2(J^^ + J'J-,B)::'^'-'' + Ja, + J'^, + 2 J,,, ; < 



(./<.2^'/;2) 

02 -T ^'^ai ^ •^12 ^ •^a2 



(104) 



2{Ja, + Ja2)\ KJ'a, + J^J; ^if "^"^^ + Ja, + J'a, + '^Ja2\B 



+ ^^n-, ^" ^a9 I ^ n 



^ a\ ^^ ^c[2 ^^ ^a2 



max 



iJo^2^J'a2) , J , jl ,r,T . ^(•■'''l '■■'»1 ) , r, 7/ 



2 J:2 + -^"1 + -^ai 



)(105) 
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2. Bulk RG rules for dynamical barriers when K = 2 
Eg 1 10 II involves a maximum over five terms 



^(ai,a„c,6i,;,.)^^^^j 



Bai ^ + Jai — Jbi + {Ja2 ^biJ'^ 

2 2 

1=1 i=l 



and Eq. 11031 involves the minimum over four terms 



(106) 



r,{ai,a2,c,bi,b2)_ Kj(a2,ai ,c,bi ,62) . K?(ai ,02 ,c,62,bl) . f2{a2,ai,c,b2,bi) 
"_R '"R '"fi '"_R 



(107) 



3. Numerical results obtained via the pool method 





60 
In L 



-5-4-3-2-10 1 2 3 4 5 



FIG. 6: Glauber dynamics of the random ferromagnetic model defined on the diamond hierarchical lattice of branching 
ratio K — 2 corresponding to the fractal dimension df = 2 : (a) Log-log plot of the averaged- value B^^^{n) and of the width 
B^idthi^) of the probability distribution Pn{B) of bulk dynamical barriers at generation n corresponding to the length L„ — 2" 



: the slope yields the dynamical exponent tp : 
distribution g{u) of Eq. 11161 



: 0.5 (b) The distribution of the rescaled barrier u : 



' B-B' 



'"W 



th(") 



) is the Gaussian 



The pool method is very useful to study renormalization rules for disordered models defined on trees [32|, |4J-|46j 
and on hierarchical lattices [^, [20, l22|, |43| • The idea of the pool method is the following : at each generation, one keeps 
the same number Mpooi of random variables to represent probability distributions. Within our present framework, 
the joint probability distribution Pp{B, J, J') of the dynamical barrier B and of the two renormalized couplings of a 
renormaUzed spin at generation n will be represented by a pool of Mpooi = 10^ triplets {Bi, Ji, J^). To construct a 
new triplet {Br, Jr, J'r) of generation (p + 1), one draws {2K) triplets {Bi, Ji, J[) within the pool of generation p and 
apply the rule of Eqs 11061 and 11071 At the last RG step p — n, one draws instead K triplets {Bi, Ji, J/) within the 
pool of generation (n — 1) and apply the rule of Eqs 11041 and 11051 to obtain the final barrier Sjf"* between the two 
ferromagnetic groundstates of the whole sample of length L„ = 2" (eq[T]). 

For the Glauber dynamics satisfying Eq. 1961 the initial condition at generation n = reads in terms of the initial 
disorder distribution p{J) of the ferromagnetic coupling reads 



P^'^I;{B, J, J') = p{J)p{J')S{B -\J- J'\) 



(108) 
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We have chosen the box distribution of width A = 1 

p( J) = 61(1 < J < 2) (109) 

On Fig. ini we present our numerical resuhs concerning the probabihty distribution of the bulk dynamical barrier 
B at generation n < 100 

pbuik(^Q^ ^ f dJ f dJ'P^^^^{B, J, J') (110) 

As explained in previous sections, this bulk dynamical barrier B characterizes the dynamics of a domain- wall crossing 
the system after its creation. We find that both the averaged value 

Bl':f^{n)= I dBBP^''^''{B) (111) 



and the width 



i?S.W ^ (^jdB B^Pi-'\B) - (i?r'W)') ' (112) 



grow with the same power-law of the length L„ = 2"^' 



Br.'*.(") ^Lt (113) 



with the dynamical exponent (see Fig. |6](a)) 



-0^0.5 (114) 

As shown on Fig. [5] (b), the corresponding rescaled barrier 

_B-Bt'Hn) 



Tjhulk ( ^\ 



(115) 



g{u) = -^e-"^ (116) 



follows the Gaussian distribution 

a(u) — — 

. 27r 

We have also computed the probability distribution of the last barrier ij'"** as a function of the generation n. 
As explained in previous sections, this last dynamical barrier B characterizes the dynamics where a domain-wall is 
created near a boundary and then crosses the system. We find that the averaged value grows linearly 

Biffin) (xL„ (117) 

as expected from the energy cost L'^' of the creation of an interface of dimension ds = df — 1 = 1 in a. space of 
dimension d/ = 2, in agreement with the result of Eq. [351 corresponding to the pure case. The width around this 
averaged value is found to scale as the width of the bulk barrier of Eq. 11131 

Bt!lM^Lt with V=^0.5 (118) 

The result -ip — 1/2 for the dynamical exponent ■0 is in agreement with the conjecture ^ — dg/2 proposed in 
our previous work [48|. In particular, ■0 = 1/2 is clearly different from the droplet exponent 9 ~ 0.299 involved 
in the statics of the random ferromagnet on the diamond lattice [49|], which coincides with the Directed Polymer 
droplet exponent 6 dp — 0.299 [21| . since the optimization of the position of the interface in the random ferro mag net 
corresponds to the optimization of the position of a directed polymer in a random medium. We refer to Ref. [i^ for 
a detailed discussion of the physical meaning of the conjecture 'ip — ds/2 with respect to other alternative proposals. 
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V. CONCLUSION 

To characterize the stochastic single-spin- flip dynamics near zero-temperature of the pure and random ferromagnetic 
Ising model on the hierarchical diamond lattice of branching ratio K with fractal dimension df = {hi{2K))/ \n2, we 
have adapted the Real Space Renormalization procedure introduced in our previous work [3(| . 

For the pure Ising model, we have obtained that the equilibrium time behaves as 

te,{L) ^ L^e^^'^"' (119) 

where ds — df — 1 is the expected interface dimension. We have computed the prefactor exponent a as a function of 
K. 

For the random ferromagnetic Ising model, we have derived the renormalization rules for dynamical barriers 
Beq{L) = (lnteq//3) near zero temperature. For the fractal dimension df = 2 (corresponding to the branching 
ratio K = 2), we have studied numerically these renormalization rules via the pool method to obtain 

BeqiL) ^ cL + L^/'^u (120) 

where u is a 0(1) Gaussian random variable of non-zero mean. The non-random term scaling as L corresponds to the 
energy-cost of the creation of an interface of dimension ds = df — 1 as in the pure case of Eq. 11191 The dynamical 
exponent tp governing the fluctuation part characterizes the barriers for the motion of a domain- wall after its creation. 
The result tp = 1/2 is in agreement with the conjecture tp = ds/2 proposed in [4^. In particular, the dynamical 
exponent ip — 1/2 is clearly different from the droplet exponent 9 ~ 0.299 involved in the statics of the random 
ferromagnet on the same lattice 49]. 

Appendix A: Renormalization rule for an effective one-dimensional dynamics 
1. First excited quantum state 

Since the exact groundstate j'i/'o > of zero energy i^o = is exactly known to be given by Eq. 1161 it is natural to 
look for the first excited state through an amplitude A{C) 



\i;,>=Y,A{C)^y^\C> (Al) 



e 2 

A[u: 

c 
Then the eigenequation for the quantum Hamiltonian of Eq. [15] 

= (H - E,)\ib, > (A2) 

can be rewritten for the amplitude A(C) as 

[Wout{C) - El] A{C) = J2 ^(^ ^ C')A{C') (A3) 

c 

2. Explicit first non-vanishing energy Ei for an effective one-dimensional dynamics 

Let us consider an effective one-dimensional dynamics between configurations (Co; Ci, C2, .., C„, C„+i) described 
by the system fEg IA3[) for 1 < i < n 

= [W{C\ -^ a_i) + W{C, -> C,+i) - ^1] A{a) ~ W{a ^ a_i)A(C,_i) - W{a ^ C,+i)A(a+i) (A4) 

and by the two boundary equations for i = and i = n + 1 

= [W{Co ^ Ci) - El] A{Co) - W{Co ^ Ci)A(Ci) 

= [W^(C„+i ^ C„) - El] A{C^+i) - W{C,,+i ^ C„)A(C„) (A5) 
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Let us assume that the intermediate configurations Ci for i — l,2,..,n have higher classical energies U{Ci) with 
respect to the two boundary configurations Cq and C„+i. Then the groundstate |i/'o > of zero energy Eq — oi Eq. 
[TBI can be approximated at low temperature by its two leading components 

1^0 > ^ - , ,^,^,[ ,^., , (e-^^(^«' |Co > +e-l^«^"-) |C„+, >) (A6) 

Then the small energy Ei can be neglected in Eq. IA4l to become for 1 < i < n 

^(Q) -P-(Q)^(C,_i) +p+(C,)A(C,+i) (A7) 

with the notations 

W{C, ^ a_i) 



P~{C^) 



W{C,-^C,^i) + W(C,^C,+i) 



^-<^') - ^(a^a5:ggJ^a.o ^'-^-'^' '^«' 



The leading components of the first excited state at low temperature 



iV-i > ^ ^ ^i(Co)|Co > +^i(C„+i)|C„+i > (A9) 



are then fixed by orthogonality with the groundstate of Eq. IA6I 



./3 



MCo) 



g-f;7(c„+i) 






so that the amplitude A{C) satisfies the boundary conditions 



A(Co) ~ -e-^lU(C^^^)-U{Co)] 

(3— S- + 00 

A(C„+i) ^ g4[c/(c„+i)-c/(Co)] (All) 

The solution of Eq. IA7I with the boundary conditions of Eq. lAllI can be obtained by recurrence [50| and reads 

A{C.) = A(Co)f|:^ ^ A[C^,,) /--f^l (A12) 

i?o(0,n) i?„+i(l,n+l) 

in terms of the Kcstcn variables [51| 

Ro{n+l,n) =0 (A13) 

RQ{n,n) — 1 



n n 



R,{k<n-l,n) =1+ Y. n^ 

1 i=m ^ ^ ' 



m—k+ 



n n 



and 



i?o(o, n) = 1 + V 1 1 — 7TT = 1 + — r^ + - + — r~\ — ? — r\ 7T\ 



Rn+i{l,0) =0 
i?„+i(l,l) =1 

k—1 m / .\ 

i?„+i(i,fc>2) =i+^n^^ 

U n ^n l^Y-fr^-W , , P-(l) , P-(l)p-(2) p-(l)p-(2)...p_(n) 
i?„+il,n+l =1+2^11 — FT = H TTT^ m — 7^ + ••■ ^ m — T^v^ TT ^^^^ 
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The energy Ei can be now computed from Eq. IA5I at the boundary Cq (or equivalently at the other boundary 
Cn+i), and using Eas lA12 land lAin one obtains 



El ~ W{Co -^ C\) 

~ W{Co ^ Ci) 
W{Co ^ Ci) 



Rn+iil,n + l) 

W{Co ^ Ci) 



Rn+i{l,n + l) 
Using Eq. IA8I and Eq. IA141 one finaUy obtains 

1 

E'i 



A{Co) 

Ro{l,n) A{Cn+l) Rn+l{lA) 

i?o(0,n) A{Co) Rn+i{l,n) 

A{Cn+l 

A{Co) 
l + e/3[t/(c„+i)-c/(Co)] 



Rn+i{l,n + l) 




1 


-U{Co)]' 


1 


-u{Co)y 






1+En 



W{Co -^ Ci) [1 + e^[^(Cn+i)-c/(Ca)]l 



1+En 






(A15) 



(A16) 



3. Renormalized amplitude Gr for an effective one-dimensional dynamics 

The renormaUzed quantum Hamiltonian is given by the projection onto the two lowest eigenstates Eq — and Ei 

H'^ff c:iEi\i:i><iPi\ (A17) 

where the first excited state 1-01 > is given by Eq lAlOl near zero temperature So Eq. IA17I becomes 



H 



eff 



El 



fi^oo e-/3'^(Co) + e-'5'^(C'"+i) 



[g-0C/(c„+i)|Co X Co\ + e-^^(^«)|C„+i X Cn+i\ 



_e-|[r/(Co)+£/(c„,.)] (|Co X C„+i| + |C„+i X Col 
So it is of the form of Eq. [TS]with only the two configurations Cq and C„+i 



13— ^+oo 



(A18) 



(A19) 



g-4[c/(c„+i)-c/(Co]|Cp X Col +e-^[^(^«)-^(^"+^l|C„+i X C„+i| - |Co X C„+i| - |a.+i >< Co 
where the renormalized amplitude reads 



Gfl(Co, C„+i) _ ^_^t/(Co) + e-^'^(C"+i) e'^'^C^o) + e/^t^C^^+i) 
Using Eq. IA16[ the final formula for the renormalized amplitude reads 

2coshf [[/(Co)-[/(a„+i)] 



(A20) 



1 

GR{Co,Cn+l) El 

g|[(7(Co)-;7(C„+i)] 

WiCo -^ Ci) 
This formula is used to obtain Eq. [551 and Eq. [SHI of the text. 



W{C, -> q-i) 



i+En 



(A21) 
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4. Example of application 

Let us now consider the case where the transition rates W{C — > C") of the effective one-dimensional problem of 
Eqs IA4I and IA5I satisfy the detailed balance form of Eq. [TU] Then we may rewrite the products as 

^ G (Co, Ci) ^4[a(c,„)+c/(c„.+i)-c/(Co)-J7(Ci)] (-^22) 

G (Cm, Cm+l) 



SO that the renormalized amplitude of Eq. IA21I reads 



g|[;7(Co)-c/(c,.+i)] 



Gfl(Co, Cn+i) G{Co, Ci)e'^VHCi)~u(Co] 



^_^y^ G{Co,Ci) ^!^iu(c^)+u(c^+i)-u(Co)-u(Ci)] 

^-^ G [CmiCm+l) 



- p-§[U(Co) + U{Cr,+i)] 



62 



[(7(C„) + (7(C„ + i)] 



e 



^ n(n n.A (^^3) 



m=0 



G (C„i, C„i+i) 



This formula is used to obtain Eq. [75] and Eq. |84]of the text. 
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